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Abstract 

We study the heavy quarkonium within the basis light-front quantization approach. We implement the 
one-gluon exchange interaction and a confining potential inspired by light-front holography. We adopt 
the holographic light-front wavefunction (LFWF) as our basis function and solve the non-perturbative 
dynamics by diagonalizing the Hamiltonian matrix. We obtain the mass spectrum for charmonium and 
bottomonium. With the obtained LFWFs, we also compute the decay constants and the charge form 
factors for selected eigenstates. The results are compared with the experimental measurements and with 
other established methods. 
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1. Introduction 

Describing hadrons from quantum chromodynamics (QCD) remains a fundamental challenge in nuclear 
physics. Inspired by the discovery of a remarkable gauge/string duality [1], holographic QCD models, 
most notably the AdS/QCD [2], have been proposed as analytic semi-classical approximations to QCD 
(for a recent review, see Ref. [3]). In light of these phenomenological successes, as well as the recent 
progress in the ab initio nuclear structure calculations [4-7], the basis light-front quantization (BLFQ) 
[8] has been developed as a non-perturbative approach to address QCD bound-state problems from first 
principles. 

BLFQ is based on the Hamiltonian formalism in light-front dynamics (LFD, [9]) in Minkowski space. 
The central task of the Hamiltonian approach is to diagonalize the QCD Hamiltonian operator, 

P+p-\ilin)=Ml\il:n). ( 1 ) 

Here = P^ ± P^ is the longitudinal momentum and the light-front quantized Hamiltonian operator, 
respectively. The eigenvalues directly produce the invariant-mass spectrum. The eigenfunctions, known 
as the light-front wavefunctions (LFWFs), play a pivotal role in the study of the hadron structures 
in deep inelastic scattering (DIS) [10] and deeply virtually Compton scattering (DVCS) [11]. In the 
Fock space expansion, Eq. (1) becomes a relativistic quantum many-body problem and can be solved by 
constructing and diagonalizing the many-body Hamiltonian matrix (see, e.g., [12] for a review). 

The advantages of LFD are made explicit by BLFQ which can employ an arbitrary single-particle 
basis subject to completeness and orthonormality. By adopting a single-particle AdS/QCD basis, BLFQ 
naturally extends the AdS/QCD LFWFs to the multi-particle Fock sectors [8]. Furthermore, this basis 
preserves all the kinematical symmetries of the full Hamiltonian [13, 14]. Such choice is in parallel 
with the no-core shell model (NCSM) used in non-relativistic quantum many-body theory [5]. State- 
of-the-art computational tools developed in the many-body theory can be used to address the QCD 
eigenvalue problem [15]. BLFQ has been applied successfully to a range of non-perturbative problems. 
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including the electron anomalous magnetic moment [16, 17], non-linear Compton scattering [18, 19] and 
the positronium spectrum [20, 21]. In this paper, we apply the BLFQ approach to the heavy quarkonium. 

Working with the full QCD Hamiltonian is a formidable task. In practice, we truncate the Fock 
space to a finite number of particles. The leading-order truncation \qq) + \ qqg) introduces the one-gluon 
exchange which produces correct short-distance physics as well as the spin-dependent interaction needed 
for the fine and hyperfine structures. The Abelian version of this interaction was extensively used in 
the literature [20, 22-25] to calculate the QED bound-state spectrum in LFD. However, the one-gluon 
exchange itself is not sufficient to reproduce the hadron spectrum since confinement is also needed. 
Holographic QCD provides an appealing approximation to confinement. 

Heavy quarkonium is an ideal laboratory for studying non-perturbative aspects of QCD and their 
interplay with the perturbative physics [26]. Conventional theoretical tools include the non-relativistic 
potential models (NRPMs) [27, 28], non-relativistic QCD (NRQCD) [29], heavy quark effective field 
theory [30], Dyson-Schwinger Equations (DSE) [31-34], and Lattice QCD [35]. The recent discoveries 
of tetraquark [36] and pentaquark [37] states have renewed interests in the theoretical investigation of 
heavy quarkonium. Extensive data on heavy quarkonium have been produced by experimental facilities, 
such as Belle, CLEO and LHC. 

Numerous light-front phenomenologies have been developed for heavy quarkonium (see e.g. [38-44] 
and the references therein). Our approach shares some similarity with these models. Yet, there are 
also major differences. First of all, our approach employs holographic QCD (confining interaction) 
and realistic LFQCD (one-gluon exchange). Secondly and most importantly, we solve quarkonium as a 
two-body bound-state problem using a Hamiltonian method that is applicable to arbitrary many-body 
bound states, once the (effective) Hamiltonian and the basis space are specified. We exploit the fact that 
BLFQ is developed as a flexible computational platform for relativistic strong interaction many-body 
bound-state problems [8, 15], designed to deal with general Hamiltonians, realistic or phenomenological. 

Our goal in this work can be simply stated: we aim to improve the light-front holographic QCD 
results [45] by including a realistic one-gluon exchange interaction. Computationally, we intend to lay 
the foundation for the extension to higher Fock sectors. 


2. Effective Hamiltonian 


2.1. Phenomenological confinement 


Our effective Hamiltonian consists of the holographic QCD Hamiltonian and the one-gluon exchange. 
We adopt the light-front AdS/QCD soft-wall (SW) Hamiltonian for the first part [46]. This simple model 
gives a reasonable description of the hadron spectrum and structures (see Ref. [45] for a review). Its 
effective “light-cone” Hamiltonian reads. 


H,. 


= p+ p- _ p2 _ 


fcj 

a;(l — x) 


+ K'^a:(l — x)r'j_, 


( 2 ) 


where, x = / P~^ is the longitudinal momentum fraction of the quark, k± = Pq± — xP± is the 

relative transverse momentum, and r± is the transverse separation of the partons. k is the strength 
of the confining potential. Note that the “light-cone Hamiltonian” has mass squared dimension, whose 
eigenvalues are the squared invariant masses. Following Brodsky and de Teramond [46], it is convenient 
to introduce the holographic coordinate Cj. = ~ x)rj_, and its conjugate q± = fcj_/\/x(l — x) = 

In light-front holography, C_l is mapped to the fifth coordinate z of the AdS space. In these 
coordinates, Hgw is a harmonic oscillator (HO), 

Hsw = Qx + ( 3 ) 


Its eigenvalues follow the Regge trajectory = 2K^{2n + \m\ + 1). Its eigenfunctions are 2D HO 
functions in the holographic variables. 




K / 


( 4 ) 
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Here q_L = |gj_|, 0 = arg< 7 j_, and L^{z) is the associated Laguerre polynomial. We adopt these functions 
as our basis. This basis has the advantage that in the many-body sector, it allows the exact factorization 
of the center-of-mass motion in the single-particle coordinates. This is a very valuable property, because 
the boson/fermion symmetrization/anti-symmetrization in the relative coordinates quickly becomes in¬ 
tractable, as the number of identical particles increases [13, 14]. For this work, however, we do not have 
identical particles in the qq sector and we will use the relative coordinate. In future extensions, as sea 
quarks and gluons are added, it may be more advantageous to adopt single-particle coordinates. 

The soft-wall Hamiltonian Eq. (3) is designed for massless quarks, and it is inherently 2-dimensional. 
For the heavy quarkonium systems, it should be modified to incorporate the quark masses and the 
longitudinal dynamics, 

^ i? 0 ON = ql+ + ^ + + Vl{x). (5) 

X \ — X 


Here is a longitudinal confining potential. Several longitudinal confining potentials have been proposed 
[47-49]. Here we propose a new longitudinal confinement which shares features with others proposed. 


Vl{x) = - 


{mq + rriq)' 


r5a;(a;(l - x)dx), 


( 6 ) 


where dx = {d/dx)^^. This term combined with the mass term from the kinetic energy forms a Sturm- 
Liouville problem. 


-dx{x{l- x)dxXi (a^)) + ^ ( 


1 / 


I — X 


— ]xi(x) = X 


)xiix) 




Xiix), 


( 7 ) 


where a = 2mq{mq + mq)l, (3 = 2mq{mq + rriq)/ The solutions of Eq. (7) are analytically known 
in terms of the Jacobi polynomial 

Xi{x) = xi“(l - x)^!^ pI°‘'^\ 2 x - 1). (8) 

and form a complete orthogonal basis on the interval [0,1]. The corresponding eigenvalues are 

a[ = {I-\-+ \{a + (3) + 1). (Z = 0,1,2, • • •) (9) 


Comparing to other forms of longitudinal confinement, our proposal implements several attractive 
features. First, the basis functions resemble the known asymptotic parton distribution ^ a;“(l — x)^ 
with a,/3 > 0 [50]. This is our primary motivation for adopting the longitudinal confinement Eq. (6). 
Second, the basis function is also analytically known, which brings numerical convenience. Third, in 
the non-relativistic limit imii{mq,mq} ^ k, the longitudinal confinement sits on equal footing with the 
transverse confinement, where together, they form a 3D harmonic oscillator potential. 


Tr _ 'IXlqlXlq 4 2 

Won — 7 i To ^ ^ 5 


(to. 


mq 


( 10 ) 


and rotational symmetry is manifest. This non-relativistic reduction also provides us an order-of- 
magnitude estimate of the model parameters for our heavy quarkonium application. Fourth, in the 
massless limit max{TO,,TOg} <C k, the longitudinal mode stays in the ground state and the longitudinal 
wavefunction Xo(a;) = const. Thus we restore the massless model of Brodsky and de Teramond^. 


2.2. One-Gluon exchange 

Following Ref. [20] (cf. [22-25]), we introduce the one-gluon exchange in LFD. In the momentum 
space, this term reads. 


(fe(i_, a;', s', s'|Eg|fej_, a:, s, s) 



47ras 

Q 9 ^S,S,S' 


, 5 /(fc_L,a;, fc^L, 


( 11 ) 


^ Note that in our normalization convention, the LFWFs differ from Brodsky et al.’s [45] by a factor y/ajT 


x). 
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where Ss,s,s',s'{k±,x,k'j_,x') = u, 
2 


+i(v^ 


(fc^,a;')7/x“s(^-L:2;)f5(-fc_L, 1 - a;)7^U5/(-fc^, 1 - a;'), and, = 
+ l(a; - x'f + Mg, is the average 

momentum transfer. We have also introduced a gluon mass fj,g to regularize the Coulomb singularity. 
This singularity is integrable and does not carry physical significance. The gluon mass is used to improve 
the numerics and will be taken small compared to other scales in this application. In Eq. (11), we have 
taken the total momentum P±^0, P+ = 1 by virtue of the boost invariance in LED. 


This one-gluon exchange introduces a logarithmic divergence [21, 23, 25, 51]. This can be seen from 
the one-gluon exchange kernel by transverse power counting. In particular, the leading power comes 
from the spin non-flip spinor matrix elements S'tltl which contain terms proportional to k'^ 

or k'^. Such terms do not vanish in the large momentum limit^. In principle, this divergence can be 
handled by a proper renormalization [52-56]. However, such a procedure is usually rather challenging or 
even impractical. In this work, we adopt a simple counterterm technique [20, 21, 23, 24, 51] by replacing 
the spinor matrix elements and with. 




S.S.S',S‘ 


(^kj _, Xj kj _, X ) y Ss^s,s',s' {kj _, X, kj _, x ) ‘2k 


I x'(l — x') 
x(l — x) 


-2k'l. 


I x(l — x) 
x'(l — x') 


(s, s, s', s' = nu or ;nt) (12) 


In its essence, this counterterm exactly removes the troubling k‘j_ and k'^ terms in the spinor matrix 
elements hence removing the UV divergence. 


3. Basis light-front quantization 

In Sect. 2, we derived the effective light-cone Hamiltonian operator for quarkonium, which reads. 


Heff — Q± + + 


fd^{x{l - x)d^) + Vg. 


X 1 — X {niq + VriqY 
In this section, we will focus on solving the eigenvalue equation non-perturbatively. 


(13) 


(14) 


to obtain the mass spectrum and the LFWFs ilj:^j{kj_,x,s,s) = (fex, a:, s, sj'i/'m/ )• Here J, P, C and 
mj are the total spin, parity, charge parity, and the magnetic projection of the state, respectively. Our 
strategy is to construct a finite-dimensional matrix from the effective Hamiltonian operator and then 
diagonalize it numerically. To do this, we first adopt a basis expansion. 


i^i,j{k±,x,s,s) = '^S, 


'mj ,Tn-\-s-\-s mj 


(n, m, /, s, s) 4>nm/(fc_L/ a/x(1 - x), x). 


(15) 


Here = (n, to, 1, s, are the LFWFs in the BLFQ basis. As mentioned above, 

the basis functions ^nmi{(lx,x) = (l)nmiq±)xii^) are generated by the confining part of the Hamiltonian 
- the kinematic energy plus the confinement. As such, the light-front holographic QCD results serve as 
our first approximation. The matrix elements of the Hamiltonian within this basis are. 


{n',m',1', s\ s'\Hes\n,m,l, s, s) = {n',m',1', s', s'\Vg\n,m,l, s, s) 


(TOq -\- 2k‘^(2ti [to] -j- / -[- k) -f 


(to, -I- rriqY 


l{l + l) 5 

nn' ^mm' 5w6ss'Sss'. (16) 


We take advantage of the conservation of the angular momentum in the transverse plane, and choose 
a particular magnetic projection toj: to -|- s -|- s = toj. In order to carry out practical calculations. 


^ From the perturbative point of view, this should be canceled by a similar contribution in the cross ladder diagram, 
which is absent from our Fock space. 
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we truncate the infinite basis space to a finite size by restricting the quantum numbers according to 
2n + \m\ + 1 < -/Vmax, I < -/Vmax controls the range of momenta covered by the harmonic oscillator 

basis as it is related to a transverse IR regulator Am ^ VV-^max and a UV regulator Ayy 
[57], where b is the HO basis scale and we have taken b = k, unless otherwise stated. Amax controls the 
basis resolution in the longitudinal direction. 

The relativistic bound states are identihed by three quantum numbers However, P is broken 

by the finite basis truncation, because the parity transformation is dynamical in LFD. Instead, one can 
exploit the so-called mirror parity Px = Rx{t^)P that is related to P [58], 

p.hC) = ( 17 ) 

Then P can be extracted from the LFWFs by, 

{-lyp = {-l)"^yiyyn,-m,l,-s,-s)yiyn,m,l,s,y. (18) 

Similarly, C can be obtained from the LFWFs by, 

n,m,l,s,s 

Rotation is also a dynamical symmetry in LFD. As a result, J is not conserved and the mass degeneracy 
for different magnetic projections mj is lifted by the Fock-space truncation. Nevertheless, in a non- 
relativistic system, such as the heavy quarkonium, the discrepancy is small and we can still extract J by 
counting the multiplicity of the nearly-degenerate mass eigenstates. It is also instructive to assign states 
the non-relativistic quantum numbers where n, S and L are the radial quantum number, the 

spin and the orbital angular momentum, respectively. L, S are related to parity and charge conjugation 
through P — (—1)^+^,C' = (—The radial quantum number n can be deduced from the mass 
hierarchy of the spectrum. The total spin S', though not an exact quantum number, can be obtained from 
its expectation value, + !)■ are constrained by the angular momentum 

addition jL — Sj < J < L -|- S. Reconstructing these quantum numbers allows us to identify the states 
and to compare with experimental data and with other methods. 


4. Numerical Results 

In this work, we focus on charmonium and bottomonium, where the fermion masses are equal (rUg = 
uiq) and heavy (m^ > k). The Hamiltonian matrix element (Eq. (16)) involves a four-dimensional 
integral (two in the radial direction, two in the longitudinal direction). They are evaluated using Gauss 
quadratures. The number of the quadrature points Nj-ad (Mfx) is taken to be at least twice AAnax (Lmax)- 
Then the obtained matrix is diagonalized using LAPACK software [59] . We fix the bottomonium coupling 
as{Mqq) — 0.25 and obtain the charmonium coupling as{Mcc) = 0.3595 using the leading-order pQCD 
evolution of strong coupling {Nf = 5,Mcc — 3.5GeV, ~ 9.5 GeV). We then ht the parameters k, 

iTT'q (jn-c, fTib) by minimizing the root-mean-squared (r.m.s.) deviation from the experimental masses [60] 
of the states below the DD or BB threshold. The details of the model parameters are summarized in 
Table 1. In these calculations, we take the regulators = 0.02 GeV ^ k, V^ax = L„iax = 8,16,24 
(A^max = 8 is available in the supplemental materials), and the number of quadrature points V^ad = 
Wfx = 64. 

Previous BLFQ study of positronium [20, 21] shows that the continuum limit V^ax —>■ oo, Lmax —>■ oo, 
fig ^ 0 can be reached through extensive calculation and successive extrapolations. However, we shall not 
investigate the Vmax, Lmax and fig extrapolations in this paper because of the numerical efforts involved 
and also the presence of several htted parameters. Comparing the results with Vmax = Lmax = 16 
and Vmax = Lmax = 24 in Table 1, where the model parameters K,mq were fitted separately and 
turned out to be very close (;< 1% change). The r.m.s. deviations from the measured spectra are 
also comparable. We also studied the trend of the mass eigenvalues with decreasing fig in the range of 
10“^ GeV < fLg < 3 X 10“^ GeV at fixed Vmax = Lmax, and found that the mass eigenvalues are well 
converged with respect to fig. 
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Table 1: Summary of the model parameters. The coupling as and the gluon mass /ig are fixed (see the 
text). The confining strength k and the quark mass are fitted with mj = 0 using the experimental 
data below the DD or BB threshold. The r.m.s. deviations for the mj = 0 spectrum, 5Mmj=o, and the 
r.m.s. average-TOj spectrum, (5M, are computed for states below the threshold. 



CXg 

Mg (GeV) 

K (GeV) 

niq (GeV) 

(MeV) 

SM (MeV) V„ax = Amax 

cc 

0.3595 

0.02 

0.950 

1.510 

64 (8 states) 

52 (8 states) 

16 

bb 

0.25 


1.491 

4.761 

56 (14 states) 

51 (14 states) 

cc 

0.3595 

0.02 

0.938 

1.522 

65 (8 states) 

52 (8 states) 

bb 

0.25 


1.490 

4.763 

54 (14 states) 

50 (14 states) 



Figure 1: Left: the charmonium spectrum in GeV; Right: the bottomonium spectrum in GeV. The 
horizontal axes is . The known states are labeled by their PDG symbols. The unknown states are 
labeled by the non-relativistic symbols. The spread of Mmj is indicated by a box and M is shown in 
dashed lines. Model parameters for our spectra shown here are given in Table 1 (A^max = A„iax = 24). 


Figure 1 shows the charmonium and bottomonium spectrum for A^max = Tmax = 24. As mentioned, 
the mass degeneracy for mj is lifted due to the violation of the rotational symmetry by the Fock sector 
truncation. We use a box to indicate the spread of masses from different mj. The r.m.s. values 

M = + Mf_j + • • • + M^j)/(2J + 1)]2 are shown as dashed bars [51]. We compare our results 

with the experimental data from the particle data group (PDG, [60], cf. [61]). The r.m.s. deviations are 
computed for charmonium (bottomonium) below the DD {BB) threshold (see Table. 1). Comparing 
with the recent results in Ref. [43], our approach improves the charmonium and bottomonium mass 
spectra from light-front holography. 

The LFWFs can be used to calculate the transition amplitudes. Here we consider the decay constants. 
These quantities are useful for computing the decay widths and constraining the Standard Model param¬ 
eters. In the non-relativistic limit, they are proportional to the wavefunctions at the origin and, therefore, 
test the short-distance physics of the model. The decay constants for a scalar S', a pseudo-scalar P, a 
axial-vector A, and a vector V states are dehned as. 


(OlV'q^y^V'lA’b)) = {Q\i^l^tp\V{pA)) = ex{p)mvfv , (A = 0,±1) 

{0\'ijjjf^ip\S{p)) = {0\ip'yf'-f^ip\A{p,X)) = e‘l{p)mAfA, (A = 0,±1) 


respectively, where e^(p) is the spin vector for the vector boson. 


e\ip) = {et{p),e^{p),ei{p)) = 


p±- 


, A = 0 


( 0 , 


2eA P 




)> 


A = ±l. 


( 20 ) 


( 21 ) 


Here = (l,±*)/-\/2. Note that due to the charge conjugation symmetry, the decay constants of 0++ 
and 1^ vanish. In LFD, the decay constants can be computed from the “-I-” component of the current 
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Figure 2: The decay constants of scalar and vector quarkonia as compared with PDG data [60] as well 
as Lattice QCD (Lattice) [62-65] and Dyson-Schwinger (DSE) [32] approaches. We extrapolate fVmax 
using second-order polynomials in and adopt the difference between the extrapolated value (Wnax 

extr.) and the fYiuax = 24 value as the uncertainty (not including systematic errors). Note that T{nD) 
are referring to the vector bottomonia T{n^Di). 


matrix elements, which read [10], 


/P,.=2v^/ 

Jo 2y'X(l — X) 

hv 


/o 2^/x{l - x) 





( 22 ) 


The decay constants from this work are plotted in Fig. 2. We also list the PDG data [60], as well as 
results from Lattice QCD (Lattice, [62-65]), and the Dyson-Schwinger equation (DSE, [32]). The PDG 
data are extracted from the dilepton decay widths Fee (for vectors) and diphoton decay widths F.^.^ (for 
pseudo-scalars). We obtain results for three successive sets of basis regulators, Nmax = f^max = 8,16, 24, 
where the parameters k and vriq are fitted to the mass spectrum separately, and as,^g are kept fixed. 
We then extrapolate A^max using simply polynomials in and estimate the uncertainty associated 

with iVniax from the difference between the extrapolated and the A^max = 24 results. While the resultant 
masses are close as we mentioned, the decay constants show noticeable residual regulator dependence. 
This may not be a surprise as the decay constant probes the short-distance physics, whereas the basis is 
chosen to emulate confinement. We expect slower convergence for the decay constants than the masses. 

From Fig. 2, our calculated decay constants are in reasonable agreement with the known experimen¬ 
tal measurements as well as Lattice and DSE results. This is encouraging since the measured decay 
constants are not used in our fits. However, compared to Lattice and DSE results, most of our results 
are systematically larger than the PDG data. This is likely due to the systematic errors of our model 
and can be improved by more realistic models for the LF Hamiltonian and/or by including higher Fock 
spaces. Fig. 2 also includes decay constants for the D-wave states. In the non-relativistic limit, these 
quantities should vanish. The small but non-vanishing D-wave decay constants in our results indicate 
the mixing of the S'-wave component, as expected from the relativistic treatment. 

The LFWFs also provide direct access to other hadronic observables. Here we study the fictitious 
charge form factor with the photon coupling only to the quark (but not the anti-quark). Of course, 
this quantity is not a physical observable. Nevertheless, it provides important insight to the system. In 
particular, this charge form factor at small momentum transfer yields the r.m.s. radius of the hadron. 


(rl) = -6^Go(02) 


Q-s-O. 


(23) 


In LFD within the impulse approximation, the form factors can be obtained from the Drell-Yan-West 
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Figure 3: The charge form factors G'o(Q^) for rjc, J/ip^ XcOi Vb, and Xbo as extrapolated from = 
Lmax = 8,16, 24 using second-order polynomials in N~l^. The difference between the extrapolated and 
the A^max = 24 values are used to quantify the uncertainty, which does not include systematic errors. 
Note that, for the scales we are showing, the charge form factor of ryt and T are on top of each other. 


Table 2: The mean squared radii (in fm^) of charmonia and bottomonia as extrapolated from iVmax = 
Lmax = 8,16, 24 using second-order polynomials in IV^ax- The difference between the extrapolated and 
the Nmax = 24 values are presented as the uncertainty, which does not include systematic errors. 


(fm^) 

«) 


/j,2 \ 

\ Xcq / 


/^2 \ 

\ r]b> 

(rl) 

/^2 \ 

\ XbO' 

(r2,) 

this work 

0.038(5) 

0.0441(8) 

0.06(1) 

0.1488(5) 

0.0146(8) 

0.0149(5) 

0.0331(8) 

0.0510(8) 

Lattice [70] 

0.063(1) 

0.066(2) 

0.095(6) 






DSE [71] 

0.048(4) 

0.052(3) 








formula within the Drell-Yan frame P''^ = P'^ [66], 


^mj ,mj 


(n)/(2P+) 


-E 


dx 


d^k± 


'o 2xil-x)J 


(24) 


where q = P' — P, and = —q^ = q\. For (pseudo) scalars, it directly produces the charge form factor 
Go(Q^) = For (axial) vector mesons, due to the violation of the rotational symmetry, there 

exists some ambiguity on finding the physical (Sachs) form factors Gq, Gi, and G 2 from [67]. 

We adopt the prescription of Grach and Kondratyuk [67] , which has been shown to be free of zero-mode 
contributions in some analytical models [68, 69]. The charge form factor according to this prescription 
reads. 


6^0 — g [(3 — 2r])Ii i + 2y/^Ii Q + _i], (25) 

where 77 = /{dM^), M is the mass of the hadron. Fig. 3 shows the charge form factors of rjc, J/i’, 

XcO, Vb, T, and Xf,o- Table 2 lists the r.m.s. radii of the first few states. From our results, the radius of 
J/ij} is close to, but slightly larger than, that of rjc- This is consistent with the Lattice and DSE results 
and can be understood from the non-relativistic point of view. Bottomonia are in general smaller than 
charmonia, a result, again, that can be drawn from the non-relativistic argument. The comparison of 
the first few available results shows that radii from our approach are in qualitative agreement with those 
of other approaches [70, 71], though ours are systematically smaller than the Lattice and DSE results. 
This observation is consistent with the trend of the decay constants. 
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5. Summary and Outlook 


We studied heavy quarkonium based on holographic QCD and a realistic one-gluon exchange on the 
light front. We proposed a longitudinal confining potential to incorporate quark mass and longitudinal 
dynamics. We solved the bound-state problem in the Basis Light-Front Quantization approach. We cal¬ 
culated the spectroscopy, decay constants and charge form factors for the charmonium and bottomonium. 
The comparison with the experimental data and results from Lattice QCD and Dyson-Schwinger equa¬ 
tion shows reasonable agreement among the available observables. As such, we improve the light-front 
holography from the first approximation to QCD. 

The LFWFs may readily be used in the study of, e.g., radiative transitions [26] and diffractive vector 
meson production in ultra-peripheral heavy ion collisions [72]. These observables will also serve as 
stringent tests of our model in different regimes. 

This model can also be applied to heavy-light and light-light systems, where the light-front holography 
provides a good initial approximation [45]. The phenomenological confinement can be improved by a 
better understanding of the string/gauge duality as well as a more complete derivation of the inter¬ 
quark potentials from various first-principle approaches to QCD [73]. Ultimately, the phenomenological 
confining interaction should be replaced by the QCD Hamiltonian. 

This work can be extended to higher Fock sectors to incorporate sea-quark and gluon degrees of 
freedom as well. The treatment of the many-body dynamics is essential for obtaining realistic predictions 
for states above the thresholds. One of the major challenges of the light-front Hamiltonian approach 
however is that explicitly including many gluons, as well as many quark-antiquark pairs, in the Fock 
space expansion quickly becomes numerically expensive. The ever increasing computational capacity and 
the progress in ab initio many-body calculations represent growing opportunities for understanding the 
strong interaction in the Basis Light-Front Quantization approach. For systems involving light quarks, 
however, the naive Fock sector truncation may not suffice, as addressing dynamical chiral symmetry 
breaking is essential [74]. The coherent basis (see, e.g., [75]) and the light-front coupled-cluster method 
[76] are two promising approaches for dealing with collective modes and dynamical symmetry breaking 
in the light-front Hamiltonian formalism. 
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Table 1: Summary of the model parameters. The coupling and the gluon mass fig are fixed. We choose the bottomonium 
Q;s(9.5GeV) = 0.25 and then obtain the charmonium coupling as(3.5GeV) ~ 0.36 from pQCD evolution of the strong 
coupling. The confining strength n and the quark mass rUg (niq = rriq) are fitted with mj = 0 using the experimental data 
below the DD or BB threshold. The r.m.s. deviations for the mj = 0 spectrum, 6Mmj=o, and the r.m.s. average-mj 
spectrum, 6M, are computed for states below the threshold. The numbers for k, rUg, SMjnj=o and 5M are rounded. 


No. 


as 

fig (GeV) 

K (GeV) 

TTlq (GeV) 

SM^g=o (MeV) 

6M (MeV) 

^ max — -^max 

I 

cc 

hi 

0.3595 

0.2500 

0.02 

0.963 

1.492 

1.492 

4.758 

65 (8 states) 

59 (14 states) 

56 (8 states) 

55 (14 states) 

8 

II 

cc 

bb 

0.3595 

0.2500 

0.02 

0.950 

1.491 

1.510 

4.761 

64 (8 states) 

56 (14 states) 

52 (8 states) 

51 (14 states) 

16 

III 

cc 

bb 

0.3595 

0.2500 

0.02 

0.938 

1.490 

1.522 

4.763 

65 (8 states) 

54 (14 states) 

52 (8 states) 

50 (14 states) 

24 


Table 2: The charmonium spectrum and decay constants (iV^ax = -^max = 8, k = 0.963 GeV,mc = 1.492 GeV). Mp^g 
quotes charmonium masses from Particle Data Group^ (PDG 2014). Mmj=o collects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass M = + M‘l_j + • • • + M^j)/(2J + 1)] ^. M^in = niin{Mmj}, M^ax = max{Mmj}. 



n^^+^Lj 


Mpdg (GeV) 

M^g=oiGeV) 

M(GeV) 

M,ai„(GeV) 

M,„ax(GeV) 

/ (GeV) 

^c(l^)* 

1% 

0-+ 

2.9836 

3.07979 

3.07979 

3.07979 

3.07979 

0.36081 

jmsr 

l^S'i 

1— 

3.096916 

3.14431 

3.10975 

3.09233 

3.14431 

0.35066 

xcoiipy 

l3Po 

0++ 

3.41475 

3.35285 

3.35285 

3.35285 

3.35285 

0. 

xciiipy 

l^Pi 

1++ 

3.51066 

3.40689 

3.44322 

3.40689 

3.46125 

0.06155 

hciipy 

VPi 

1+- 

3.52538 

3.49972 

3.49546 

3.49333 

3.49972 

0. 

xc2iipr 

1^P2 

2++ 

3.5562 

3.51397 

3.50420 

3.48788 

3.51557 


rfc{2S)* 

2% 

0-+ 

3.6394 

3.69728 

3.69728 

3.69728 

3.69728 

0.29471 

ip{2S)* 

2^Si 

1— 

3.686109 

3.72362 

3.68679 

3.66824 

3.72362 

0.30062 

b(3770) 

l^Di 

1— 

3.77315 

3.74542 

3.72814 

3.71946 

3.74542 

0.07773 


I^D2 

2— 


3.75804 

3.77571 

3.75804 

3.78275 



VD2 

2-+ 


3.81879 

3.79676 

3.78678 

3.81879 



l^Da 

3— 


3.83039 

3.80230 

3.77322 

3.83039 


Xcoi^P) 

2 3Po 

0++ 

3.9184 

3.88304 

3.88304 

3.88304 

3.88304 

0. 

V(3872) 

2^Pi 

1++ 

3.87169 

3.91125 

3.95026 

3.91125 

3.96962 

0.07521 

V(3900)± 

2^Pi 

1+- 

3.8887 

4.01287 

4.00042 

3.99419 

4.01287 

0. 

Xc2(2F) 

2^P2 

2++ 

3.9272 

4.01768 

3.99424 

3.96949 

4.01768 



3% 

0"+ 


4.20176 

4.20176 

4.20176 

4.20176 

0.25915 

b(4040) 


1— 

4.039 

4.21367 

4.17769 

4.15958 

4.21367 

0.26696 


* States whose masses are used to fit the parameters k and rric in this work. 

^ K.A. Olive et al. (Particle Data Group), Chin. Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 


Table 3: The charmoninm spectrum and decay constants (A^max = Tmax = 16, k = 0.950 GeV, me = 1.510 GeV). Mpdg 
quotes charmonium masses from Particle Data Group^ (PDG 2014). Mmj=o collects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass M = [{M^ j + M^_ j H— • + M^j)/(2J + 1)] ^. M^in = niin{Mmj}, M^ax = max{Mmj}. 


_ n^^+^Lj jPg Mpdg(GeV) M(GeV) M^i^jGeV) M^^^jGeV) /(GeV) 

r;e(lS')* 1% 0-+ 2.9836 3.07306 3.07306 3.07306 3.07306 0.40055 

J/b(15')* l^S'i 1— 3.096916 3.16284 3.12261 3.10230 3.16284 0.37249 



Xco(lP)* 

l^Po 

0++ 

3.41475 

3.33761 

3.33761 

3.33761 

3.33761 

0. 

Xcl(lP)* 

I^PI 

1++ 

3.51066 

3.40861 

3.44986 

3.40861 

3.47030 

0.05899 

hciipy 

I^PI 

1+- 

3.52538 

3.51013 

3.50708 

3.50555 

3.51013 

0. 

Xc2{IPY 

1^P2 

2++ 

3.5562 

3.52568 

3.51702 

3.50142 

3.52823 


ric{2S)* 

2% 

0"+ 

3.6394 

3.67956 

3.67956 

3.67956 

3.67956 

0.32243 

y{2S)* 

2 3S'i 

1— 

3.686109 

3.72234 

3.67897 

3.65710 

3.72234 

0.33323 

V'(3770) 


1— 

3.77315 

3.74999 

3.73280 

3.72418 

3.74999 

0.08694 


I^D2 

2— 


3.76292 

3.78028 

3.76292 

3.78602 



1^D2 

2"+ 


3.81867 

3.80099 

3.79281 

3.81867 



13^3 

3— 


3.83233 

3.80674 

3.77982 

3.83233 


Xco(2F) 

2^Po 

0++ 

3.9184 

3.84589 

3.84589 

3.84589 

3.84589 

0. 

X(3872) 

2^Pi 

1++ 

3.87169 

3.89320 

3.93699 

3.89320 

3.95870 

0.09212 

X(3900)± 

2^Pi 

1+- 

3.8887 

4.00585 

3.99503 

3.98960 

4.00585 

0. 

Xc2{2P) 

2^P2 

2++ 

3.9272 

4.01273 

3.99120 

3.96669 

4.01273 



3% 

0-+ 


4.17128 

4.17128 

4.17128 

4.17128 

0.29910 

V'(4040) 

3 3S'i 

1~" 

4.039 

4.19738 

4.15054 

4.12692 

4.19738 

0.28728 

* States whose masses are used to fit the parameters k, and nic 

in this work. 



^ K.A. Olive et al. (Particle Data Group), Chin. 

Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 


Table 4: The charmonium spectrum and decay constants (A^max 

— -^max ~ 

= 24,k = 0.938 GeV, me = 1.522 GeV). Mpdg 

quotes charmonium 

masses from Particle Data Group^ (PDG 2014). Mmj=o collects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass 

M = [{M^j + M^_j + 

••• + M2^)/(2J+1)]A 

^min = min{M„j}, M^ax = 

max{M„j}. 

71 



Mpdg(GeV) M„^=o(GeV) 

M(GeV) 

M^i„(GeV) 

M„,ax(GeV) 

/ (GeV) 

iciisy 

1% 

0-+ 

2.9836 

3.06877 

3.06877 

3.06877 

3.06877 

0.42302 

Jiy{isy 

l^S'i 

1— 

3.096916 

3.17612 

3.13355 

3.11205 

3.17612 

0.37965 

xcYipy 

l^Po 

0++ 

3.41475 

3.32860 

3.32860 

3.32860 

3.32860 

0. 

xcYipy 

l^Pi 

1++ 

3.51066 

3.41041 

3.45460 

3.41041 

3.47649 

0.05406 

hyipy 

l^Pi 

1+- 

3.52538 

3.51612 

3.51375 

3.51256 

3.51612 

0. 

xc2{ipy 

1^P2 

2++ 

3.5562 

3.53182 

3.52396 

3.50908 

3.53486 


ric{2Sy 

2% 

0-+ 

3.6394 

3.66753 

3.66753 

3.66753 

3.66753 

0.33044 

y{2sy 

2 3S'i 

1— 

3.686109 

3.72126 

3.67578 

3.65283 

3.72126 

0.34403 

V'(3770) 


1— 

3.77315 

3.75121 

3.73492 

3.72675 

3.75121 

0.08943 


13^2 

2— 


3.76390 

3.78070 

3.76390 

3.78569 



1^D2 

2"+ 


3.81578 

3.80081 

3.79382 

3.81578 



13D3 

3— 


3.83024 

3.80650 

3.78128 

3.83024 


Xco{2P) 

2 3Po 

0++ 

3.9184 

3.82321 

3.82321 

3.82321 

3.82321 

0. 

A:(3872) 

2 3Pi 

1++ 

3.87169 

3.88279 

3.92842 

3.88279 

3.95105 

0.09723 

A:(3900)± 

2 3Pi 

1+- 

3.8887 

3.99890 

3.98910 

3.98419 

3.99890 

0. 

Xc2(2P) 

2 3^2 

2++ 

3.9272 

4.00641 

3.98650 

3.96288 

4.00641 



3% 

0-+ 


4.15002 

4.15002 

4.15002 

4.15002 

0.30702 

V'(4040) 

3 3S'i 

1— 

4.039 

4.18663 

4.16370 

4.15218 

4.18663 

0.28008 

* States whose masses are used to fit the parameters k and me 

in this work. 



^ K.A. Olive et al. (Particle Data Group), Chin. 

Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 


Table 5: The bottomonium spectrum and decay constants (A^ma> 

: — -^max 

= 8,k = 1.492 GeV, mb = 4.758 GeV). Mpdg 

quotes bottomonium 

. masses 

from Particle Data Group^ (PDG 2014). Mmj= 

0 coliects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass 

M = [{M^j + M^_j + 

••• + M2^)/(2J+1)]A 

^min = min{M„j}, M^ax = 

max{M„j}. 

n 

25+1^^ 

JP^ 

Mpdg(GeV) M^,=o(GeV) 

M(GeV) 

M^i„(GeV) 

M^ax(GeV) 

/(GeV) 

vbYsy 


0-+ 

9.398 

9.50231 

9.50231 

9.50231 

9.50231 

0.46979 

T{isy 

l3S'i 

1— 

9.4603 

9.52295 

9.50957 

9.50288 

9.52295 

0.45714 

xboiipy 

l"Po 

0++ 

9.85944 

9.82548 

9.82548 

9.82548 

9.82548 

0. 

hbilpy 

l^Pl 

1+- 

9.8993 

9.86487 

9.86623 

9.86487 

9.86691 

0. 

xbiiipy 

l^Pl 

1++ 

9.89278 

9.83976 

9.84875 

9.83976 

9.85325 

0.02550 

Xb2iipy 

13P2 

2++ 

9.91221 

9.86946 

9.86803 

9.86272 

9.87262 


ribi^Sy 

21^0 

0-+ 

9.999 

10.04033 

10.04033 

10.04033 

10.04033 

0.39332 

T{2sy 

2 3S'i 

1— 

10.02326 

10.05169 

10.03727 

10.03005 

10.05169 

0.39680 




l^Di 

1— 


10.11386 

10.11010 

10.10822 

10.11386 

0.01723 

T(1D)* 

I^D2 

2— 

10.1637 

10.11785 

10.12319 

10.11785 

10.12463 



I^D2 

2-+ 


10.13580 

10.13294 

10.13177 

10.13580 



13^3 

3 — 


10.13542 

10.13272 

10.12568 

10.13597 


xboi^py 

2 3Po 

0++ 

10.2325 

10.28144 

10.28144 

10.28144 

10.28144 

0. 

hb{2P)* 

2^Pi 

1+- 

10.2598 

10.32862 

10.32696 

10.32614 

10.32862 

0. 

Xbi{2Pr 

2^Pi 

1++ 

10.25546 

10.29200 

10.30551 

10.29200 

10.31226 

0.03913 

Xb2{‘2P)* 

2^P2 

2++ 

10.26865 

10.33087 

10.32467 

10.31497 

10.33127 



I^F2 

2++ 


10.36309 

10.35457 

10.34957 

10.36309 



13^3 

3++ 


10.36561 

10.36731 

10.36419 

10.37097 



11^3 

3+- 


10.38030 

10.37436 

10.37182 

10.38030 



13^4 

4++ 


10.38223 

10.37372 

10.36290 

10.38223 



3% 

0-+ 


10.51184 

10.51184 

10.51184 

10.51184 

0.35815 

T(35')* 

3 3S'i 

1— 

10.3552 

10.51945 

10.50394 

10.49617 

10.51945 

0.36402 


2^Di 

1— 


10.54918 

10.54435 

10.54193 

10.54918 

0.02425 


2 3^2 

2— 


10.55200 

10.56022 

10.55200 

10.56478 



2^D2 

2-+ 


10.58138 

10.57365 

10.57153 

10.58138 



2 3^3 

3 — 


10.58271 

10.57122 

10.55865 

10.58271 



1^03 

3 — 


10.59956 

10.58520 

10.57695 

10.59956 



13G4 

4 — 


10.60157 

10.59853 

10.58970 

10.60684 



1^4 

4“+ 


10.61546 

10.60405 

10.59818 

10.61546 



1^5 

5 — 


10.61732 

10.60322 

10.58792 

10.61732 


X&?(3P) 

3 3Po 

0++ 

10.5794 

10.72027 

10.72027 

10.72027 

10.72027 

0. 


3 3Pi 

1++ 


10.72703 

10.74293 

10.72703 

10.75087 

0.04392 


3 3Pi 

1+- 


10.77096 

10.76538 

10.76259 

10.77096 

0. 


3 3P2 

2++ 


10.77131 

10.75943 

10.74662 

10.77131 



4% 

0”+ 


10.95999 

10.95999 

10.95999 

10.95999 

0.31175 

T(4S') 

4 3S'i 

1— 

10.5794 

10.96411 

10.95056 

10.94377 

10.96411 

0.29659 

* States whose 

masses are used to fit the parameters k and mb 

in this work. 



^ K.A. Olive et al. (Particle Data Group), Chin. 

Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 


Table 6: The bottomonium spectrum and decay constants (Amax 

— -^max — 

16, K = 1.491 GeV,mb = 4.761 GeV). Mpdg 

quotes bottomonium masses 

from Particle Data Group^ (PDG 2014). Mmj=o 

collects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass 

M = [ 

'{Mlj + M^_j + 

• • • + MIj)/{2J + 1)] G = min{M„,}, = 

max{Mmj}. 


n^^+^Lj 


Mpdg(GeV) M^,=o(GeV) 

M(GeV) 

M^i„(GeV) 

M^ax(GeV) 

/(GeV) 


1% 

0-+ 

9.398 

9.49520 

9.49520 

9.49520 

9.49520 

0.53938 

T(IS')* 

l3S'i 

1— 

9.4603 

9.52317 

9.50657 

9.49826 

9.52317 

0.50967 

xboiipy 

l3Po 

0++ 

9.85944 

9.82577 

9.82577 

9.82577 

9.82577 

0. 

hyipy 

I'Pi 

1+- 

9.8993 

9.86971 

9.87106 

9.86971 

9.87173 

0. 

xbiiipy 

13P3 

1++ 

9.89278 

9.84213 

9.85228 

9.84213 

9.85735 

0.01536 

Xb2{ipy 

13P2 

2++ 

9.91221 

9.87471 

9.87324 

9.86777 

9.87798 


rib{2Sy 

2% 

0-+ 

9.999 

10.03214 

10.03214 

10.03214 

10.03214 

0.46047 

T{ 2 sy 

2 3S'i 

1— 

10.02326 

10.04925 

10.03120 

10.02216 

10.04925 

0.45700 


13^1 

1— 


10.11883 

10.11464 

10.11254 

10.11883 

0.01827 

TilD)* 

13^2 

2— 

10.1637 

10.12300 

10.12834 

10.12300 

10.12990 



11^2 

2-+ 


10.14099 

10.13816 

10.13710 

10.14099 



13^3 

3 — 


10.14346 

10.13838 

10.13081 

10.14346 


xboi2py 

2 3Po 

0++ 

10.2325 

10.27462 

10.27462 

10.27462 

10.27462 

0. 

hb{2Py 

2 3Pi 

1+- 

10.2598 

10.33064 

10.32905 

10.32825 

10.33064 

0. 

Xbii2py 

2 3Pi 

1++ 

10.25546 

10.28884 

10.30469 

10.28884 

10.31260 

0.03509 

Xb2{‘2py 

2 3P2 

2++ 

10.26865 

10.33363 

10.32740 

10.31715 

10.33452 



13^2 

2++ 


10.36803 

10.35942 

10.35435 

10.36803 



13^3 

3++ 


10.37062 

10.37225 

10.36921 

10.37571 




3+- 


10.38509 

10.37936 

10.37683 

10.38509 



13^4 

4++ 


10.38710 

10.37873 

10.36791 

10.38710 



31^0 

0-+ 


10.49843 

10.49843 

10.49843 

10.49843 

0.43906 

T{isy 

335 '! 

1— 

10.3552 

10.51176 

10.49099 

10.48059 

10.51176 

0.44228 


2^Di 

1— 


10.55174 

10.54460 

10.54102 

10.55174 

0.02241 


2 3^2 

2— 


10.55509 

10.56334 

10.55509 

10.56747 



2^02 

2-+ 


10.58458 

10.57689 

10.57415 

10.58458 



2 3^3 

3 — 


10.58627 

10.57435 

10.56074 

10.58627 





l^Gs 

3 — 


10.60407 

10.58968 

10.58107 

10.60407 



l^Gd 

4 — 


10.60613 

10.58410 

10.54711 

10.61116 



l/b'd 

4“+ 


10.61971 

10.60863 

10.60269 

10.61971 



l^Gs 

5 — 


10.62171 

10.60776 

10.59234 

10.62171 


Xb?{3P) 


0 ++ 

10.5794 

10.70393 

10.70393 

10.70393 

10.70393 

0 . 


3^Pi 

1 ++ 


10.71532 

10.73571 

10.71532 

10.74589 

0.04951 


3iPi 

1 +- 


10.76921 

10.76330 

10.76034 

10.76921 

0 . 


3^P2 

2 ++ 


10.77023 

10.75765 

10.74316 

10.77023 



41^0 

0 "+ 


10.93561 

10.93561 

10.93561 

10.93561 

0.41974 

T(4S') 

435 '! 

1 — 

10.5794 

10.94581 

10.92267 

10.91109 

10.94581 

0.42446 


* States whose masses are used to fit the parameters k and mb in this work. 

^ K.A. Olive et al. (Particle Data Group), Chin. Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 

Table 7: The bottomonium spectrum and decay constants (IVmax = Amax = 24, k = 1.490 GeV,m{, = 4.763 GeV). Mpdg 
quotes bottomonium masses from Particle Data Group^ (PDG 2014). M^j^o collects the masses from the mj = 0 sector. 

The r.m.s. average-mj mass M = [(M^j + M^_j H— • + M^j)/(2J + 1)] ^. Mmin = niin{Mmj}, Mmax = max{Mmj}. 


n^s+^Lj Mpdg(GeV) M™,=o(GeV) M(GeV) M^i„(GeV) M^ax(GeV) /(GeV) 



l^So 

0-+ 

9.398 

9.49182 

9.49182 

9.49182 

9.49182 

0.57499 

T(IS')* 

l^Si 

1— 

9.4603 

9.52401 

9.50588 

9.49680 

9.52401 

0.53223 

XboilP)* 

l^Po 

0++ 

9.85944 

9.82570 

9.82570 

9.82570 

9.82570 

0. 

XbiilPy 

l^Pi 

1++ 

9.89278 

9.84297 

9.85371 

9.84297 

9.85907 

0.00766 

hbiipy 

I'Pi 

1+- 

9.8993 

9.87168 

9.87301 

9.87168 

9.87367 

0. 

Xb2iipr 

13P2 

2++ 

9.91221 

9.87680 

9.87533 

9.86982 

9.88009 


ribi^sy 

21^0 

0-+ 

9.999 

10.02867 

10.02867 

10.02867 

10.02867 

0.49181 

T{2S)* 

2 3S'i 

1— 

10.02326 

10.04909 

10.02950 

10.01969 

10.04909 

0.48152 


l^Di 

1— 


10.12060 

10.11636 

10.11424 

10.12060 

0.01891 

T{1D)* 

13^2 

2— 

10.1637 

10.12479 

10.13013 

10.12479 

10.13172 



13^2 

2-+ 


10.14277 

10.13996 

10.13891 

10.14277 



l^Ds 

3 — 


10.14525 

10.14019 

10.13264 

10.14525 


Xboi^P)* 

2 3Po 

0++ 

10.2325 

10.27171 

10.27171 

10.27171 

10.27171 

0. 

Xbii^Py 

2^Pi 

1++ 

10.25546 

10.28761 

10.30453 

10.28761 

10.31298 

0.02967 

hb{2P) 

2^Pi 

1+- 

10.2598 

10.33154 

10.32993 

10.32912 

10.33154 

0. 

Xb2{2Py 

2 3P2 

2++ 

10.26865 

10.33473 

10.32853 

10.31821 

10.33575 



13^2 

2++ 


10.36949 

10.36092 

10.35586 

10.36949 



13^3 

3++ 


10.37208 

10.37371 

10.37070 

10.37714 



l^Fs 

3+- 


10.38651 

10.38082 

10.37831 

10.38651 



13^4 

4 ++ 


10.38852 

10.38020 

10.36941 

10.38852 



33^0 

0 -+ 


10.49364 

10.49364 

10.49364 

10.49364 

0.47294 

T{3sy 

335 '! 

1 — 

10.3552 

10.51028 

10.48765 

10.47632 

10.51028 

0.47155 


2^Di 

1 — 


10.55258 

10.54524 

10.54157 

10.55258 

0.02352 


2 3^2 

2— 


10.55602 

10.56429 

10.55602 

10.56834 



2 3^2 

2-+ 


10.58552 

10.57790 

10.57517 

10.58552 



2 3^3 

3 — 


10.58730 

10.57542 

10.56180 

10.58730 



13G3 

3 — 


10.60517 

10.59085 

10.58226 

10.60517 



13G4 

4 — 


10.60722 

10.60423 

10.59558 

10.61220 



13G4 

4“+ 


10.62070 

10.60975 

10.60385 

10.62070 



l^Gs 

5 — 


10.63009 

10.60955 

10.59352 

10.63009 


Xb?{SP) 

3 3Po 

0 ++ 

10.5794 

10.69830 

10.69830 

10.69830 

10.69830 

0 . 


3 3Pi 

1 ++ 


10.71193 

10.73392 

10.71193 

10.74490 

0.04780 


3 3Pi 

1 +- 


10.76906 

10.76307 

10.76007 

10.76906 

0 . 


3 3P2 

2 ++ 


10.77028 

10.75775 

10.74307 

10.77028 



43^0 

0 ”+ 


10.92944 

10.92944 

10.92944 

10.92944 

0.45811 

T{4S) 

435 '! 

1 — 

10.5794 

10.94290 

10.91748 

10.90475 

10.94290 

0.45977 


* States whose masses are used to fit the parameters k, and mb in this work. 

^ K.A. Olive et al. (Particle Data Group), Chin. Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 



Table 8: The masses, decay constants for pseudo-scalar and (axial-)vector quarkonia^. See Table 1 for model parameters. The 
masses are taken from the r.m.s. average-mj values M. The decay constants are computed using the light-front wavefunction 
of the mj = 0 sector. The decay constants are also extrapolated (extr.) using simple polynomials in iV“ax- 


mass M (MeV), iV^ax = 8,16, 24 decay constant (MeV), iVmax = 8,16, 24 




PDG2 

1 (8) 

II (16) 

III (24) 

PDG2»b 

I (8) 

II (16) 

III (24) 

extr. 


Vc(lS) 

2983.6 

3080 

3073 

3069 

330±13 

361 

401 

423 

482 


Ve(2S) 

3639.4 

3697 

3680 

3668 

000 +72 
^oo_i04 

295 

322 

330 

345 

0 -+ 

Vb(lS) 

9398.0 

9502 

9495 

9492 


470 

539 

575 

665 


r/b(2S) 

9999.0 

10040 

10032 

10029 


393 

460 

492 

568 


Vb(3S) 


10512 

10498 

10494 


358 

439 

473 

551 


Ji^ 

3096.916 

3110 

3123 

3134 

407±5 

351 

372 

380 

394 


^(2S) 

3686.109 

3687 

3679 

3676 

290±2 

301 

333 

344 

366 


ii(3770) 

3773.155 

3728 

3733 

3735 

97.7±3 

78 

87 

89 

94 

1 — 

T(1S) 

9460.3 

9510 

9507 

9506 

689±5 

457 

510 

532 

585 

T(2S) 

10023.3 

10037 

10031 

10030 

479±4 

397 

457 

482 

537 


T(3S) 

10355.2 

10504 

10491 

10488 

414±4 

364 

442 

472 

535 




10110 

10115 

10116 


17 

18 

19 

21 


T(2^Di) 


10544 

10545 

10545 


24 

22 

24 

28 


Xcl(lP) 

3510.66 

3443 

3450 

3455 


62 

59 

54 

38 

1 ++ 

XbiilP) 

9892.78 

9849 

9852 

9854 


26 

15 

8 

0 * 

Xbi(2P) 

10255.46 

10306 

10305 

10305 


39 

35 

30 

13 


Xbi(3P) 


10743 

10736 

10734 


44 

50 

48 

39 


* The naive extrapolation yields a negative value (—14 MeV), which is corrected to 0, as the quantity should be positive. 

^ The decay constants of O"*"’*" (scalar meson) and 1"' {h axial vector meson) are exactly zero due to charge conjugation 
symmetry. 

^ K.A. Olive et al. (Particle Data Group), Chin. Phys. C 38, 090001 (2014); [http://pdg.lbl.gov]. 


“ For pseudo scalars (0 +), the decay constants are extracted from the diphoton decay width according to^ 

TT 2 

>77 I-^P7 (^) I > 

4 nip 

where eg is the charge number of the constituent quark (eg = 2/3 for charm quark and eg = —1/3 for bottom 
quark), mp is the mass of the pseudo scalar, Ofem = cXemiQ) is the QED running coupling. For charmonium, we 
take aem(3.5 GeV) ~ 1/134; for bottomonium, we take aem(9.5 GeV) ~ 1/132. The higher order pQCD correction 
is assumed to be included in the non-perturbative dynamics of the bound states. 

^ For vector mesons (1 ), the decay constants are extracted from the dilepton decay width according to^“^, 


47r 


- V—>-e+e- — 


2 2 
ena^„ 


‘my 


where my is the mass of the vector meson. 

3 S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985). 

^ C.T.H. Davies, C. McNeile,E. Eollana, G.P. Lepage, H. Na, and J. Shigemitsu (HPQCD collaboration), Phys. Rev. D 
82, 114504 (2010). 



